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Abstract 

I consider the expansion of transcendental functions in a small parameter around rational 
numbers. This includes in particular the expansion around half-integer values. I present 
algorithms which are suitable for an implementation within a symbolic computer algebra 
system. The method is an extension of the technique of nested sums. The algorithms allow in 
addition the evaluation of binomial sums, inverse binomial sums and generalizations thereof. 
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1 Introduction 



The expansion of higher transcendental functions [1,2] occurs frequently in many areas of sci- 
ence. In particular, one encounters these functions in the calculation of higher order corrections 
to scattering processes in particle physics. In a previous publication we considered the expan- 
sion of transcendental functions in a small parameter around integer values [3,4]. The restriction 
to integer values is in general sufficient for the evaluation of loop integrals arising in massless 
quantum field theories. However, the inclusion of particle masses in loop integrals [5-13] or 
the evaluation of phase space integrals [14-20] can lead to half-integer values. It is therefore 
desirable to extend the algorithm of [3,4] to include at least half-integer values. Here I report on 
algorithms for the expansion of transcendental functions around rational number p/q, where p 
and q are integers. In particular this includes the half-integer case. 

Each term in the expansion is expressed through multiple polylogarithms [21,22]. Compared 
to the pure integer case, the extension to rational numbers p/q introduces naturally the ^-th roots 
of unity in the arguments of the polylogarithms. All algorithms are based on manipulations of a 
special form of nested sums. These nested sums are generalizations of Euler-Zagier sums [23,24] 
or harmonic sums [25-29]. The algorithms presented here can be implemented into a symbolic 
computer algebra system like Form [27,30] or GiNaC [31]. 

As a spin-off, the methods presented here allow the evaluation of binomial sums [5,6], inverse 
binomial sums [5,7-9,32,33] and generalizations thereof. Inverse binomial sums are sometimes 
evaluated with the help of log-sine integrals [10, 11]. In an appendix I compare the log-sine 
approach with the one presented here. 

This paper is organized as follows: Sect. |2lrecalls the definition and main properties of nested 
sums and multiple polylogarithms. It is a brief summary of ref. [3]. Sect. |3]introduces roots of 
unity and gives the basic algorithms for the expansion around rational numbers. Sect. |4| treats 
binomial sums and generalizations thereof. Sect. |5] deals with inverse binomial sums and gen- 
eralizations thereof. Sect. I^gives some simple applications to massive loop integrals and phase 
space integrals. Finally, sect. |7]contains a summary and the conclusions. An appendix compares 
this approach to log-sine integrals and collects some important relations for polylogarithms of 
low weight. 

2 A summary of known properties of nested sums 

In this section I shortly review properties of particular forms of nested sums, which are called 
Z-sums and 5-sums. Details can be found in [3]. Z-sums are defined by 

Z{n;mi,...,mk\xu...,Xk) = £ T^IT'-'T^ 
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and form a Hopf algebra. If the sums go to Infinity (n = oo) the Z-sums are identical to multiple 
poly logarithms^ [21]: 

Z{oo;mi,...,mk;Xi,...,Xk) = Umi,...,mk{xh ■■■,Xk). (2) 

For XI = ... = = 1 the definition reduces to the Euler-Zagier sums [23, 24]: 

Z(R;mi,...,mfe;l,...,l) = Zmj,...,mjn). (3) 

For n = oo and xi = ... = x;^ = 1 the sum is a multiple ^- value [22]: 

Z(oo;mi,...,mA;; !,...,!) = ^mu-,mk- (4) 

The multiple polylogarithms contain as the notation already suggests as subsets the classical 
polylogarithms Li„(x) [34], as well as Nielsen's generalized polylogarithms [35] 

Sn,p{^) = Li„+i,i,...^i(x,l^^^^, (5) 

p-l 

the harmonic polylogarithms [36] 

Hm\,...,mii{x) = Li»ii,...,>Wifc('^; J; 1 ) (6) 

k-1 

and two-dimensional harmonic polylogarithms [37]. The usefulness of the Z-sums lies in the 
fact, that they interpolate between multiple polylogarithms and Euler-Zagier sums. In addition, 
the interpolation is compatible with the algebra structure. 

In addition to Z-sums, it is sometimes useful to introduce as well iS-sums. S-sums are defined 

by 

5(n;mi,...,m/,;xi,...,x*,) = £ -ji^.-'-jki 

n>h>i2>->ik>l ^ ^ 

and form an algebra. The S-sums reduce for xi = ... = x;^ = 1 (and positive mi) to harmonic 
sums [27]: 

5(n;mi,...,mjt;l,...,l) = (8) 

The S-sums are closely related to the Z-sums, the difference being the upper summation bound- 
ary for the nested sums: {i— 1) for Z-sums, i for 5-sums. It is advantageous to introduce both 
Z-sums and 5-sums, since some properties are more naturally expressed in terms of Z-sums while 
others are more naturally expressed in terms of S-sums. One can easily convert from one notation 
to the other. 

^Note that we use here the reversed notation for multiple polylogs and multiple zeta values as compared to 
refs. [3,21]. 
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Basic manipulations involving nested sums are: 
Conversion: 



Multiplication: 



Z{n\...) 

S{n;...) Z(n;...). (9) 



Z(n;...) Z(n;...) ^ Z(n;...), 

S{n-...)S{n-...) ^ S{n-...). (10) 



Convolution: 



E ^Z(/-l;...)^;^Z(n-/-l;...) Z(n-1;...). (11) 



Conjugation: 



;=1 

Conjugation and convolution 



|;( • ^^O;-) - %;-)• (12) 



• J(-i)'];;;^0";--0^;;^ - 5(n-i;...). (13) 

These algorithms are described in [3]. It is worth to recall some technical steps to evaluate the 
conjugation in eq. (fT^ . since the same pattern of steps will be used in sect. HI To evaluate the 
conjugation it is convenient to introduce yet another type of sum as follows: 

N N N 

B{n;N;mi,...,mk;xi,...,Xk) = L ••• L fn^-^t- (1"^) 

These 5-sums can be used to express S-sums with upper summation limit n in terms of 5-sums 
with upper summation limit A'^: 

S{n\mi,...,mk\xi,...,Xk) = S{N\mi,...,mk\xi, ...,Xk) - S{N\m2,...,mk\X2, ...,Xk)B{n\N\mi\x\) 
+S{N\m3,...,mk;x3,...,Xk)B{n;N;mi,m2;xuX2) - ... + {-l)''B{n;N;mi,...,mk;xi,...,Xk). 

(15) 

Eq. (fT5l allows also to express a 5-sum recursively in terms of 5-sums S{N;...) and 5(n; ...): 

B{n;N;mi,...,mj,;xi,...,Xk) = {-lfS{n\mi,...mk\xi,...,Xk) 

~{-\fS{N;mi,...mk;xi,...,Xk) + {~l)''S{N\m2,...mk\X2,...,Xk)B{n;N-,mi;xi) 

- ... + {-if S{N\mk\Xk)B{n\N\mi, ...,mk-\\x\, ...,Xk-\) (16) 
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Finally, it is convenient to introduce raising and lowering operators as follows: 

1 



(x+)'"-l = ^In-W, 
^ ' ml 



ax 



(17) 



It is understood that in the second line only functions which are integrable at jc = are considered. 
With the help of raising operators a 5-sum B{n;oo; ...) may be expressed as follows: 



5(n;oo;mi,...,m^;jci,...,x^) 



(Xk j (Xk-1 j ...(Xi j - 



m\ Xk ^k—l^k 



XI... Xk 

-Xk 1 -Xk-lXk 1 -XI. ..Xk 

Some important integrals related to the raising operators are: 

^ 1 



{xi...xky\ (18) 



X+[1-(1-XiX2)"] = £t l-{l-XiX2y 



!=1 



XiX2 
1 -XiX2 



^ 1 



[l-(l-JCoXiJC2)"] =-(l-Xo)" J^T ( - 

7=1 ' 



1 - (1 -^0^1X2) 



+(i-(i-^or)i;^+x 
7=1 ' 



1 1 (-^1^2)^(1 - (1 --yo)") 



1 — X1X2 



(19) 



To evaluate the conjugation in eq. dT^ one first converts the 5-sum to 5-sums and introduces 
then the raising operators. This allows to perform all sums explicitly and one is left with the 
integrals corresponding to the raising operators. With the help of eq. (fT9t one can systematically 
perform these integrals and convert them back into nested sums. 

The basic manipulations in eqs. Q to (fT3l) are the building blocks to reduce the following 
two generic types of sums to single Z-sums. 
Type A: 



^ x' r{i + ai+bie) r{i + ak + bke)^,. ^ . 
~[{i + crr{i + ci+di£) r{i + Ck + dk£) 



(20) 



Type B: 



n-l 



jc' r{i + ai+bie) r{i + ak + bk£) 



y - 

^tt (« + c)'"r(z + ci+Ji£) r{i + ck + dke) 



Z(z + o- l,mi,...,m/,.JCi,...,JC/) 



r{n - i + a\ + b\£) T{n - i + a'^, + b',^,t) 



{n-i + c'y r(n - i + c\ +d[z)"' T{n - i + 4, + 4,8) 
yiZ[n — i + o — l,mj, ...,m;,,Xi, ...,X;/J 



(21) 
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Here, all aj, a'j, cj and c'j are integers, c, c', are nonnegative integers and o, o' are integers. For 
sums of type A the upper summation limit n may extend to Infinity. In the rest of the paper I relax 
the condition on a^, a'y, Cj and c'j to allow for rational numbers. Ref. [3] contains in addition two 
generic types (labelled type C and D in [3]) involving a conjugation. These types do not allow an 
extension to rational numbers along the lines of this paper. The conjugation is however related to 
two other important types of sums, which are generalizations of binomial and inverse binomial 
sums. They are treated in sect. HJand sect. |51of this paper. 



3 Algorithms for the expansion around rational numbers 

In this section I extend the algorithms for the expansion of transcendental functions towards the 
expansion around rational numbers. Subsection B. ll introduces roots of unity. Roots of unity are 
useful for the refinement algorithm in subsection l3.21 which allows to express a 5-sum ...) as 
a combination of 5-sums, whose upper summation limit is an integer multiple of n. Subsection 
I3.3l treats the expansion of Euler's Gamma function around rational numbers. Finally, in subsec- 
tion l3.4l all pieces are assembled and the algorithms for the expansion of functions of type A and 
B around rational numbers are given. A restriction on these algorithms is given by the fact, that 
rational numbers have to appear in the same place in the numerator and in the denominator. This 
restriction is relaxed in sect. |4|(a rational number only in the numerator) and sect. |5l(a rational 
number only in the denominator). 



3.1 Roots of unity 

We define a short-hand notation for the roots of unity: 

= exp(?=i^). (22) 

Here, (^r is a positive integer and p is a non-negative integer. We will need a few properties of the 
^-th roots of unity. Powers of the ^-th roots of unity are periodic modulo q: 

{rlY^' = {rPy- (23) 
Sums of powers of the ^-th roots of unity yield: 

^ ~ 10, m^Omod^. ^ ' 

If m = mod q, the proof of the relation is trivial. In the case m^O mod q we may assume that 
< m <q (due to the periodicity). Then ^ 1 and we have 

E«r=E W = -r^=<'- (25) 

p=0 p=0 ^ 'q 
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From eq. (l24b we obtain immediately 



^ h I 0' otherwise. ^ 

Since this sum occurs frequently, we introduce the notation 

= -1(4) • (27) 

Eq. (l26t will be useful to convert sums with upper summation limit n to sums with upper sum- 
mation limit (q-n). Note that dp^q{m) is idempotent: 

5p,g(m) ■ 5p,^(m) = hp^q{m) for all integer m. (28) 

Some examples for the half-integer case are: 

5o,2w = ^[i+(-in, 

Klin) = ^[l-(-in. (29) 



3.2 Refinements of ^-sums 



A 5-sum 5(n;mi, ...) with upper summation limit n can be expressed as a combination of 
5-sums with upper summation limit (q-n), where qis a positive integer. The algorithm proceeds 
recursively in the depth of the 5-sum. For the empty sum one has 

S{n) = S{q-n). (30) 

For a 5- sum of the form 

5(n;mi,m2,...;xi,JC2,...) = V -p5(/;m2, ...;JC2, ...) (31) 

the algorithm converts first the subsum S(/;m2, ...;x2, •■•) to a combination of subsums S{q-i;...). 
Finally, the outermost sum is converted according to 

tisiq-i;...)=q-f^do,q{i)l-(x'/A'sii;...) = q'^^^^ (32) 

i=l ' ;=1 ^ ^ ^ p=07=l ^ 

This completes the algorithm for the conversion of 5-sums S{n; ...) to S{q ■ n;...). As an example 
we have 

5(n;l, l;xi,X2) = 5(2n; 1, 1; ^^) +5(2n; 1, 1; -v/^) 

-|-5(2n; 1, 1;-^^, V^) +5'(2n; 1, 1; -^/xi)- (33) 

Note that as a consequence of the refinement algorithm one obtains relations like 

U,n{x^) = 2"'-i[Li,„(x)+U„(-x)]. (34) 
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3.3 Expansion of the Gamma function 

For the expansion of the Gamma function around positive integer values one has the well-known 
formula 

- r(„^,)e.p(-|e'(^..w). (35) 

For the expansion around rational numbers one finds 



1 J \ 1 

r(n + \-£-] ( 1^-1 



^^^exp ( -1 £ £ Z^^-^S{q.n-Xr\) . (36) 

Here, n and q are positive integers and p is an integer with < p < q and gcd(p,^) = 1. For the 
case p = and q = I this reduces to the formula eq. (B5t . Eq. is derived as follows: One 
starts from the expansion of the logarithm of the Gamma function: 

lnr(^n + l-^ + £^ = lnr(^n + l-^^ + f;^V'"i)(^n + l-^^. (37) 

Here \\f^^^^'> (x) is the polygamma function defined as 

^\f^'-'\x) = ^lnr(x). (38) 
From the recurrence formula for the polygamma function 

\\r^'-'\x+\)-\\f^'-'>ix) = -r{k){--] (39) 



one obtains 



'(x+i)-V'"i)(x) = -r(^)(-^)' 



The sum in the last term of eq. (HUl) is then rewritten in terms of 5-sums with upper summation 
limit (q-n) with the help of eq. (l26l) . 
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3.4 Expansion around rational numbers 

In this section we generalize the algorithms A and B, which expand the transcendental sums in 
eq. (l20b and eq. (1211) around integers to algorithms for the expansion around rational numbers. 
In the following we will always assume that qj is a positive integer and that pj is an integer with 
< Pj < qj and gcd(pj,^j) = 1. We restrict ourselves here to ratios of Gamma function of the 
form 



r(n + c,-g+J,£) 



(41) 



where n, Uj and bj are integers. Here the same fraction Pj/qj occurs in the numerator and the 
denominator. In products of these ratios, different fractions are allowed, like for example in 

r(n + ai - ^ + &ie) r(n + a2 - i + ^2e) 

— = — (42) 

r(n + ci-i+Ji£) Y{n + C2-l + d2Z)' 

The restriction to ratios of the form ensures that the prefactors Y{n +1 — p/q) from the 
expansion of the Gamma function in eq. cancel between the numerator and the denominator. 
The restriction in eq. (HTT) will be relaxed in sect. Eland sect. |5] We now consider sums of the 
type A and B, where we allow the substitutions 

_Pl _Pi 
Qj Qj 

r' r' -EI 

CI j ^ CI j . , C j ^ C j . , 

Qj Qj 

PO I I p'q rAi\ 
C C , C C j-^ (43) 

?0 Qq 

in eqs. (l20b - (|2l])- For example, a sum of type A involving rational numbers is of the form 

where cq is a non-negative integer, o and all ay, Cj are integers. When dealing with rational 
numbers it is more convenient to work with 5-sums instead of Z-sums and we replaced the Z- 
sum in eq. (EUb by a 5-sum in eq. (EH). Due to the conversion algorithm the two formulations 
are equivalent. 

The algorithms for the expansion in £ of these sums starts by reducing the offsets o and o' in 
the subsums 5(/ + o;mi, ...;;ci, ...) and S{n — i + o';m[, ...■,x[^...) to zero. Then, using the identity 

r{x+i) = xr{x) (45) 
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for the Gamma function, the ratios of the Gamma functions are brought to the form 



— % or . (46) 

They are then expanded in £, using eq. (l36b . This yields 5-sums with upper summation limit 
(qj ■ i) or {q'j ■ {n — i)). Now, let q = \cm{qQ,qi, ...,qk) be the least common multiple of qo, 
q\, qk and let q' = lcm(^Q,^j, ...,^^) be the least common multiple of ^q, q[, q'^ Using 
the refinement algorithm we can convert any occuring 5-sum S{qji;...) to 5-sums with upper 
summation limit (q-i). Similar any 5-sum S{q'j{n — /);...) is converted to 5-sums S{q'{n — ?');...). 
Products of 5-sums are then converted into single 5-sums with the help of the multiplication 
algorithm. After partial fractioning one arrives at the following forms 



Type A: £- -7^5(^/;...) 



V 90 



Type B: £ ^S{qi-...)S {q\n - i);...) , (47) 



i=i (i + c-^ 

where qo divides q. For sums of type A one then reduces the offset c to zero and one arrives at 
sums of the form 

ItT^ow^^-), (48) 



Using eq. (l26t this sum can be written as 



n i qn \x 

i=l qo> i=\ ' 



= I Ux'l'^) i:i^^5(/-f;.;...), (49) 

where p = poq/qo is an integer with < p < q. Finally, reducing the offset p of the subsum 
S{i + p;...) to zero, one arrives at 5-sums with upper summation limit (qn). 

For sums of type B one first refines the subsums 5(^/; ...) and 5 {q'{n — /);...) to S{qi;...) and 
S{q{n — /);...), respectively. Here q = \cm{q,q') is the least common multiple of q and q'. The 
next step consists in rewriting 

n— 1 

L 7 ^S{qi;...)S{q{n-i);...) = 



^"V/^ £ ^yn Sii + p;...)Siqn-i-p;...), (50) 



qn—q 



Alq' 



i=\ 



{i + qc) 
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which brings us back to the integer case. 

This completes the necessary modifications for the extension towards the expansion around 
rational numbers for the algorithm A and B. 



4 Binomial sums and generalizations 



Here we study sums of the form 



1 



r(l-E],tl r(n + l) n'"o 



5(n;mi,...,m^;xi,...,Xyt). 



(51) 



This relaxes the condition (HIT) and allows one unbalanced fraction 

r(n+i- 



r(n + i) ■ 

A special case of the form in eq. (BTT) are binomial sums 



«=i 



2n 



S{n\mi,...,mk\xi,...,Xk), 



which are obtained by setting p = \, q = 2, xq = Az'm (BTT) and by noting the identity 



In 
n 



(52) 



(53) 



(54) 



r(i)r(n + i)- 

Binomial sums have been studied in [5, 6]. The evaluation of the sums in eq. (BTT) follows the 
same pattern of steps as the evaluation of the conjugation in eq. dT^ . as given in ref. [3]. Due 
to eq. (fTSt we may replace in eq. (BTT) the 5-sum S{n;...) with 5-sums 5(oo;...) and 5-sums 
5(n;oo; ...) and it is sufficient to study sums of the form 



rU + 1 



Y(\-P^ti r(n + l) n-o 



B{n;oo;mi,...,mk;xi,...,Xk) 



mo -^fc ^k—l^k 



-(xk+) (xk-i+) ...(xi+) (xo+) Y3^Y3 



X\...Xk 



-Xk-\Xk \-Xl...Xk 



X 



— 1 — ^ 
\-{\-XQXl...Xk) ^ 



(55) 



To derive the r.h.s of eq. ( 1551 we first introduced the raising operators x+, then performed 
explicitly all geometric sums from 5(n;oo;0, ...,0;xi, ...,Xk) and finally performed the remaining 
sum with the help of the hypergeometric summation formula 



r(i-f) -1 



T{n + l) 



l-(l-jc) V t 



(56) 



11 



We then perform succesivly the integrations corresponding to the raising operators. Due to the 
appearance of rational numbers in the exponent, the set in eq. ( fT^ is no longer sufficient and has 
to be supplemented with additional equations. Nevertheless, the principle stays the same: Each 
integration preserves the structure, such that multiple integrations can be performed iteratively. 
To reduce mo we have 

n=l ^ 



l-(l-xoxi) ^ ^' 



^O^q-p,q{n) l-{l-X0Xl)l =q^q-p,q{n)Y,-[^0,q{i)-^q-p,q{^)] 1 " ( 1 " .^0-^1 ) ' 



i=l 



A 1 



+q8q-p,q{n) -5q-p,q{i) 1 - (1 -XoXl)'' 



i=l 



(57) 



To reduce mi, nij, nii^we have in the general case xq 1 

, XiX2 



- 1-- 

-ll l-(l-XoXiX2) ^ " 

1 -XiX2 

oo 

E ^ [hq{n)-^q-pM] ( Y 



-^(l-^o) ' 



X 



n=l 



1 



■xo 



1 -(1 ~XoXiX2y' 



+ 



l-(l-Xo) V ' 



N 



y - {xiXj^+xt ^^^^ {XIX2) 

n ^1 — jcijc? 



A' 



Xj8q-p_q{n)- 



X\X2 



-X\X2 L 



,^1" ^ XXX2 

1 - (1 -XoJClJC2)« = -^(1 -^o)*5^-p,<7(« 



l-(l-;co) y ^ 



A 1 



1 



■Xq 



1 - (1 -X0XiJC2)« 



^ 1 ^1X2 r «" 

^'^q-p,q{n) \-{\-XqY Y^-{xxX2)' ^hq-p^q{n)x\— {xxX2f l-(l-Xo)* . 

z=l ^ L 



(58) 



In the special case xo = 1 we have instead: 



X\X2 
1 -X\X2 



\-{\-X,X2) ^ 



1 



P 



1-(1-Xi.^2) ^ 



(59) 



-X\X2 



.^1^2 
1 -.x:iX2 



1 -(1 -X1X2)'' 



N 



1 -(1 -jci.x:2)'? + ^ T {xiX2y+x^ {X1X2) 



TV 
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Eqs. (l58t and (l59l) introduce an arbitrary integer A'^. It is worth to note that the eqs. (l58t and (l59l) 
hold for any integer A^. We will take the limit 0° in the end. It is clear that in this limit terms 
of the form 



can be neglected. Some examples are: 



r^+^mx^ ^ y ^ (60) 



4n E ^7^- = 21n2 - 21n f 1 + v^' 

wTY E , -2 = 2 ln21n 1 + VT^) +ln21n 1 - VT^) - Liu - v^T^, 1 
-Lin (v^r^,-l)-Li2(-l)-(ln2)2] , 



rG)„tij^ir(n+i) n J ' y \-x j ' ' 'V i-^ 

-2 Liu (^VT^,--jl=^ -2Un (^-v^r^, -^=L=^ +4Li2(-/r^) 

-4Li2(-l)+21nl±^lniiVp. (61) 

V 1— 

5 Inverse binomial sums and generalizations 

In this section I consider sums with one unbalanced rational number in the denominator. Exam- 
ples of these type are sums of the form 

11=1 r ( Hi " 



^/ «,=1 rin, 

Wi-l "2 n/t-2-l v-"*^-! 



X V -^2 V (-62^ 



Note the slightly different structure of — n/) in the denominators of the subsums. Only in 
the innermost sum a factor appears in the denominator. Changing the summation variables 
according to 

ni= ji+j2 + --- + jk, n2 = j2 + --- + jh «fc = jk, (63) 
this sum is equivalent to 

r ( 1 - - ) L - L -r. : — ; — ^7^~r^- — — ■ ^^^^ 
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A special case of this form are inverse binomial sums 



E 



1 

n 



(65) 



which are obtained by setting p = I, q = 2, a = b = I, k = 1 and xi = z/4 in eq. (l64l) . Inverse 
binomial sums have been studied in [5, 7-9, 32, 33] and have been related in the literature to 
log-sine integrals [10, 11]. Here I follow a different path for the evaluation of inverse binomial 
sums. The method presented here is closely connected to the transformation 



—X 



(66) 



whose inverse transformation is again given by 



—X 



To start with the evaluation of sums of type (l64b one replaces all factors 1 /j by 

j X-^oTU+l+l) 



(67) 



(68) 



The original sum is recovered as the ?i — > limit of the regularized sum. The introduction of the 
regularization in X allows us to extend the lower summation boundary of all sums from 1 to 0. 
This can be done recursively according to 



I/(7) = -/(0) + I/(7)- 
/(O) corresponds to a sum of lower depth. We therefore deal with sums of the type 



/ \ oo oo 



7r 



r(ji + ...+A+<3) 
1 jwr(ji + ...+A+Z7-| 

.i.(.>«)-...(.......r[rg^]'". 

For these sums we have the following integral representation: 



rUk + i+i) 



rrik 



q 



dun! 



a-l 



{\-u) 1 



dtxt[ 



dt:, t 



r(z>-a- 

{\-Uti...ti^Xiy^ {\-Uti^ + l...ti^XlX2y^ ... {\-Uti^_^ + i...ti^Xl...Xk) 



-1 



(69) 



(70) 



(71) 
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where 

/i=mi, i2 = mi+m2, ik = nii + ... +mk. 

For the integration variable u we now perform the substitution 

M ^ 1 — M, 

and we obtain for eq. (1711) 



(72) 



(73) 



r(i-f)r(a) ^ £ r(ji + ...+A+z>-a -f) | 
r(^&-a-|) 71=0 j,=o r(^ji + ...+jk+b- 



dtit[ 



dtij 



Ik ip. 



q 



{-tl...ti^Xiy^ {l-tl...ti^Xl) ^1 ^ ... [-tii^_^ + l...tii^Xl...Xky'' (l -?/j_j + l...?,-^JCl...JC;t) 

r(l-f)r(a) ^ "^ir(ni+Z,-a-f 



£\ 1 



A,-l 



dt- 1^^^ 







-iik-i+rik-l 



{-t\...ti^X\)"^ "2 (1 -ti...ti^X\] 

Note that the substitution eq. (TTSt induces the transformation x — x/(l — x) as follows: 



[1-(1-m);c1 



'l-x) 



I 



—X 



(74) 



(75) 



As a sideremark it is worth to note that as special case one obtains by this procedure the trans- 
formation of the hypergeoemetric function 



2Fi{a,b;c;x) = (l-x) " 2F1 i a,c - b;c; 



\-x 



(1 -x) 2F\ ( c-a,b;c 



l-x 



The purpose of the substitution in eq. (17^ is to change the ratio of Gamma functions from 
r{ji + ... + jk + a) r[j, + ...+jk + b-a-(^ 



r(ji + ...+A+Z7-f 



to 



r(ji + ...+A+z,-f 



(76) 



In the last form, the rational number p/q appears both in the numerator and the denominator 
and can therefore be treated with the methods discussed in sect. |3l It remains to perform the 
integration over ti, ...,ti^. These integration are recursively done according to the formula 



j dtitly^ {-tQtixY {l-tQtixY 



r(n + a) f. r(fc + c) ^ ^ . ^\-k-c 
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As a net result we obtain a rooted tree with side-branches, which can be expanded in X and 
converted to nested sums with the help of the algorithms discussed in ref. [3] and in sect. EI In 
the final result all poles in X, cancel and one can extract the ^I'^-term, which yields the evaluation 
of the inverse binomial sum eq. (l62b . Some examples are: 



r 
r 
r 



2j,hr{n + ^)n ^ 1+x' 

I J = -Liii (X, 1) +LiH (X, -1) -Liu (-%, 1) +Liii -1) , 



-Liu K 1) -Liu (-"0, 1) +Liu (\), -1) +Liu (-"U, -L 



-1 



Liu ( ) +Liu 1 %, 



^l-Li„ 

X 



-X,^)-Liu 



-%ln(l -xy)\n 



1+%' 



(77) 



where % = ^J—x/{l —x) and \) = ^J—xy/{l —xy). 



6 Applications 

In this section I give some applications of the techniques described in the previous sections 
relevant to the calculation of Feynman loop integrals with massive particles or to phase space 
integrals. The first example concerns a one-loop triangle with a uniform internal mass mi=m2 = 
m^ = m and two vanishing external momenta p\ = p\ = 0. The third external momentum = 
is kept arbitrary. Such diagrams are for example relevant to Higgs physics. Specific examples 
are Higgs production via gg ^ H or Higgs decay into two photons // — > 7y. In both cases the 
interaction proceeds via an internal top quark loop. It is well known that the corresponding 
integral in dimensional regularization is proportional to [7, 38] 

Here £ = (4 — D)/2 denotes the deviation from four space-time dimensions. The expansion of 
the hypergeometric function to sufficient high order in £ is a non-trivial task [7,9]. The expansion 
can be accomplished systematically with the methods presented here. First one notices that there 
is one unbalanced half-integer number in the denominator. With the results of sect. |5]one can 
show that 

3F2 (^l+aie, 1 -fa2£, 1 +a3£; ^+b£,C + ce; = {l-zy^^"'" 

r(§+te)r(c+c£) ^r{n + i+ai£) f -z Y 



r(l-fai£)r(l-Fa3£)r(i-F(Z7-a2)£)r(C-l + (c-a3)£)„to T{n + C + c£) \l-z, 



16 



" r(7 + l+a3e)r(j + i + (Z7-a2)e) r{n- j + C-l + {c-a3)£) 



r{j + l+be) 



r(n-j + i; 



(79) 



This yields a nested sum where all half-integer numbers are balanced in the numerator and in 
the denominator. This expression can now be expanded systematically in £ with the algorithms 
described in sect. El For the example mentioned above one finds 



3F2( 1,1, l+£; -,2;z 



1 



Liii 
+2e 



,1 +Liii 



Li 



11 



-z 



Li 



111 



-Li 



111 



1-z 

-z 



,1,1 +Liii] 



-,1,1 



1-z' 
-Liui 



Li 



11 



-,-1 



—z 



1,-1 



1 



1,-1 



+ o(e2) . 



(80) 



As a second example I discuss briefly phase space integrals. Analytic results for phase space 
integrals are needed for example in the calculation of higher order corrections for jet observables 
as integrals over approximation terms within the subtraction method [14, 17, 18,39]. An example 
is the integral over a subtraction term, which approximates the emission of a soft gluon from a 
heavy quark pair. This integral is given by [14] 



'^ee(^o,e) 



-2£ 



Jrr-2e-i(l-r)-^(l-ror)-i 



ds{l 



[(2(l-ror)-(l-ro))(l-5o5)-i-(l-ro)(l-5o5)-2] 



where 



'ro(l -r) 



l-ror 

This integral corresponds to the following triple sum 



'^'ee(^o,e) 



r(i 



r(£) 



-2e 



OO CO oo 



III(-iy(i+(-iy>S 

;=0 j=Ok=0 



k+j/2 



Tji + J + 1 ) r(/ + e) Tjk - 2e)r( j72 + 1 - e)r(fe + j/2) 

r(/ + J + 2 - £) r(/ +i)r{k+ \)r{j/2)r{k + j/2 + 1 - se) 

2-(l-ro)(2 + y):^ 
Using the methods of sect. |3]this sum can be expanded systematically in £. The result is 



'J^ee(^o,e) 



. ll+ro l + v^ 
2 ./rr, 1 - 



21nro — In 



2 f^ + Vni 



+ 



1 .J^ + V^o' 



:ln 



(81) 



(82) 



(83) 
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-Li2 



Li2 (v^) - Li2 (- V^) + 2 Li2 



+ ln21n 
+0(£). 



+ U2 



1 + 



1 + 

+ ln(l-v^) In 



-Li2 

1 + 



1 + 

2 



2Li2 



21nroln 



+ Li2 



+ iln2(l + y^)-lln2(l-0^) 



(84) 



7 Summary and conclusions 

In this paper I reported on algorithm which allow the expansion of certain transcendental func- 
tions in a small parameter around rational values. These algorithms extend the ones for the 
expansion around integer values and are based on the manipulation of specific forms of nested 
sums. Of particular importance is the case of the expansion around half-integer values. This 
case occurs frequently in the calculation of radiative corrections in quantum field theories with 
massive particles. The methods presented in this paper allow a systematic approach for the cal- 
culation of these integrals. 



A Log-sine integrals 

Inverse binomial sums are often expressed in terms of log-sine integrals. In this appendix I briefly 
summarize the results from the literature. The following inverse binomial sum can be evaluated 
with elementary functions as follows: 



E 



r{n+l) x" 2v^arcsin(v^) 



(85) 



This result agrees with the one given in eq. (1771 . In the literature, evaluations of inverse binomial 
sums of higher weights are given in terms of log-sine functions [5,7-11, 32, 33] : 



T{n + \) z" 



m—2 



i-iy 



where = 2arcsin ^ and the log-sine functions are defined by 



(86) 



Lsf (0) 



dQ' (0yin^' 



-k-l 



2 sin 



0' 



(87) 



By analytic continuation the log- sine functions are then related to poly logarithms [9]. A simple 
example is given by 



Lsf (0) = Cl2(0), 



(88) 
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involving the Clausen function CI2. The Clausen functions Cl„ are given in terms of poly loga- 
rithms by 

r ^.[U„{e'^)-Un{e-'^)], neven, 
Cl„(e) = <^ (89) 
[ i[Li„(e'«)+Li„(e-'9)], nodd. 

B Relations for polylogarithms 

Multiple polylogarithms of low weight can be expressed in terms of ordinary logarithms and 
polylogarithms. Relations relevant to the examples in this paper are 

Lii (;c) = — ln(l— ;c), 
Liii(x,l) = ^\n^{l-x), 

Liii(jc,-1) = ^^2-^ln^2-ln(l-x)ln(l+Jc)+ln21n(l+Jc)-Li2 



Un{x,y) = ln(l-x)ln(l-3;)+Li2f-^VLi2^~^^^"''^ 



1 , 

Liiii(;£:, 1,1) = --ln-'(l-x), 



Liiii(x,-1,-1) = i^2ln2-J^3-^ln32-iln(l-x)ln2(l+x) + ^ln21n2(l+x) 

Z O D Z Z 

-ln(l +x)Li2 +Li3 . (90) 

More relations can be found in [37,40]. 
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